Structure of entanglement in regulated Lorentz invariant field theories 
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Regulated Lorentz invariant quantum field theories satisfy an area law for the entanglement entropy S of a 
spatial subregion in the ground state in d > 1 spatial dimensions; nevertheless, the full density matrix contains 
many more than e non-zero eigenvalues. We ask how well the state of a subregion R in the ground state of such 
a theory can be approximated when keeping only the e"^ largest eigenvalues of the reduced density matrix of R. 
We argue that by taking the region R big enough, we can always ensure that keeping roughly e states leads to 
bounded eiTor in trace norm even for subregions in gapless ground states. We support these general arguments 
with an explicit computation of the error in a half-space geometry for a free scalar field in any dimension. 
Along the way we show that the Renyi entropy of a ball in the ground state of any conformal field theory at 
small Renyi parameter is controlled by the conventional thermal entropy density at low temperatures. We also 
reobtain and generalize some old results relevant to DMRG on the decay of Schmidt coefficients of intervals in 
one dimensional ground states. Finally, we discuss the role of the regulator, the insensitivity of our arguments 
to the precise ultraviolet physics, and the role of adiabatic continuity in our results. 









I 

o 
o 



> 

O 

o 
en 



X 



I. INTRODUCTION 

Finding useful representations of quantum states is a task 
of fundamental importance in the physics of quantum matter 
It is a problem that we must master to fully understand quan- 
tum matter, but it is also exceedingly practical: such under- 
standing is almost certain to lead to powerful techniques for 
studying the physics of concrete and relevant physical models. 
For example, the power of DMRG [ 1] to efficiently solve one 
dimensional systems arises from an efficient representation of 
ground states in one dimension known as matrix product states 
(MPS) (see Ref. [2] for a recent review). 

Entanglement plays a basic role in the search for such 
efficient representations. In concrete terms, the existence 
of an area law [3, 4] for the entanglement entropy of spa- 
tial subregions in many systems suggests a generalization 
of matrix product states called tensor network states (or 
PEPS/MERA/etc.) [5-7]. In a tensor network state the am- 
plitude for some configuration of physical degrees of freedom 
is computed by contracting a large local network of tensors. 
The locality of the network then guarantees that an area law 
for entanglement entropy is obeyed. 

To make these ideas more precise, consider a bipaitite quan- 
tum system divided into regions R and R with the whole sys- 
tem in its ground state (with some unspecified local Hamilto- 
nian). The Schmidt decomposition of the full quantum state 
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in terms of the Schmidt coefficients A^ and orthonormal 
Schmidt vectors \A, i) (A = R, R). The reduced density ma- 
trix of R is 

Pfl-tr^dV^X^D^E^'l^'*)^^'*!' (1-2) 

i 

and the entanglement entropy is 

SiR) = -tVRipiilnpR) = -^A.lnA,. (1.3) 



We say the the entanglement entropy of a region R of linear 
size L in the ground state of a local d-dimensional Hamilto- 
nian satisfies an area law if the entropy scales like S{R) ^ 
\dR\ ^ L'^^^. Since e'^ provides a rough measure of the di- 
mension of the support of pR, i.e. the number of non-zero 
Xi, it seems interesting to consider a class of states having 
the property that the density matrix of a region R has support 
bounded by e"^'^^' for all R. Tensor network states provide 
just such a class since entanglement between R and its com- 
plement R is only generated by bonds that cut the boundary 
of i?. Hence the entropy of region R is bounded by \dR\ logx 
where x is the "bond dimension" of the tensors forming the 
tensor network, i.e. the indices of the tensors run from 1 to at 
most X- 

However, we hasten to emphasize that the above discussion 
is heuristic at best since the physical density matrix pR will 
not in general have only e'^ non-zero eigenvalues. Indeed, we 
expect to have many more than e'^ non-zero eigenvalues of 
PR, SO a truncation of the support of pR entails some error 
One precise sense in which e'^ is the right number of eigen- 
values is in the asymptotic setting where one can show that 



® can be well approximated by a state with support on 

'Ms 



States provided M is large enough [8]. 

In this work we are interested in the following question: 
what is the error in directly approximating pR (not p^") by a 
state with only ~ e'^ non-zero eigenvalues. We are interested 
in this question in part because it has bearing on the suitabil- 
ity of tensor network states to represent ground states of lo- 
cal Hamiltonians. If, as we show, e'^ states suffice, then ten- 
sor network states, even with constant bond dimension, have 
enough degrees of freedom to represent ground states of local 
Hamiltonians very precisely (to small constant error in trace 
norm). Such a result represents non-trivial evidence favoring 
the tensor network variational class beyond evidence obtained 
from numerical simulations of particular models. 

Our questions are analogous to those asked in Refs. [9, 10] 
except that we focus on ground states of local Hamiltonians 
instead of generic (Haar random) states. Since we cannot an- 
swer our question in general we focus on an interesting class 
of systems that can be called regulated Lorentz invariant quan- 



turn field theories. These are lattice models or otherwise reg- 
ulated models with a finite number of degrees of freedom that 
approximate Lorentz invariant field theories at low energy. A 
review of a subclass of such models relevant for the physics 
of strong interactions can be found in Ref. [11, 12]. More 
generally, examples of these models abound in nature, includ- 
ing electrons at half-filling in graphene (provided we neglect 
or screen Coulomb interactions) and quantum magnets in the 
Ising universality class at their critical point. All known topo- 
logical phases in two dimensions also have representatives in 
this class [13] (although we certainly do not claim that all 
topological phases are Lorentz invariant), so our tools apply 
there as well. 

All the models we consider (in d > 1 spatial dimensions) 
satisfy an area law for the entanglement entropy, S{R) ~ 
\dR\, even if they are gapless. We use powerful geometri- 
cal tools that give access to the entire spectrum of the reduced 
density pR for special subregions R (balls and half-spaces) 
to compute errors [14-16]. These tools have a long history 
in the physics of quantum fields in curved spacetime, black 
holes, and quantum gravity (see Ref. [17] for a nice intro- 
duction), and now we are putting them to work describing the 
structure of entanglement in certain kinds of quantum matter 
The half-space results apply to any Lorentz invariant theory 
while the ball results require the additional presence of con- 
formal symmetry. Our notion of error is usually provided by 
the trace norm ||A||t = tr{VA^A) which has a natural statis- 
tical interpretation. Requiring that two states be close in trace 
norm, e.g. ||pi — /92||t < e, is a stringent condition (see Ap- 
pendix A for a discussion), but we will see that it is possible 
with rather modest resources. 

We now summarize our results. We show that the Renyi 
entropy Sa, defined for general region R by 



S^iR) = 
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■log(tr(pS)), 



(1.4) 



in the special case where i? is a d-ball in the ground state 
of any conformal field theory (CFT) obeys an area law when 
a < 1 (but not not too small), and that the ground state Renyi 
entropy as a —> is controlled by the thermal entropy den- 
sity of the CFT. This is a perhaps surprising result, namely 
that ground state properties may be used to determine low 
temperature thermal properties. Using a result of Ref. [18] 
this immediately implies that bounded error can be achieved 
keeping only e^° states. In detail, we can achieve error e in 
trace norm by keeping e'^"^'' states such that log e < — i^^^(5 
(0 < a < 1). Note that this bound does not directly apply 
to the entanglement entropy since a = 1 is excluded, but we 
argue that a similar bound holds even when a = 1 for the 
models we consider We show that that same is true for half- 
space regions in any Lorentz invariant theory, and we analyze 
in detail the case of a free scalar field showing explicitly that 
the error can be made small for a half-space. Another product 
of our analysis is a bound on the decay of Schmidt coeffi- 
cients even in gapless states. As we discuss in more detail be- 
low, this enables us to recover and generalize old results about 
the decay of Schmidt coefficients in the context on DMRG on 
gapped spin chains. 



Let us elaborate briefly on the issue of tensor network 
states. We prove that the density matrices in question can be 
well approximated by their e'^ largest eigenvalues, but we do 
not prove that the full system has an efficient tensor network 
representation. One way to rephrase this issue is to say that we 
consider the density matrix in the basis in which it is diagonal 
(thus effectively rendering the problem classical), but this ba- 
sis may not be easily constructible from a tensor network state 
by tracing out some regions. Our result is necessary but not 
sufficient if we want a constant (in system size) bond dimen- 
sion X- 

The remainder of this work is organized as follows. First, 
we briefly review existing results related to our question. Sec- 
ond, we discuss our tools and give some general arguments. 
Third, we verify our general arguments with detailed error c al- 
culations for the simplest field theory: the free scalar We 
conclude with some broader discussion of our arguments and 
possible extensions. In an appendix we discuss different met- 
rics for comparing quantum states and record a pathological 
example illustrating the need, in principle, for our considera- 
tions in the main text. 



II. REVIEW OF EXISTING RESULTS 

Before turning to our main arguments, we will briefly re- 
view some of the most relevant existing results regarding our 
main questions. We always denote the density matrix of re- 
gion R by pb. and the entropy by S{R), and we define the 
"entanglement Hamiltonian" Kf{ via pfj = cxp {—Kp;). The 
Schmidt coefficients defined above are then related to "entan- 



glement energies" Ei (spectrum of Kji) via A^ 



The 



set of all the Ei is called the "entanglement spectrum" [29]. 
The area law for entanglement entropy says that the entropy 
S{R) ^ \dR\, a result that naturally leads us to suspect that 
Kpi may define an effective local Hamiltonian in one lower 
dimension. We will see to what extent this is true below. We 
will speak about efficiency loosely below, but roughly what 
we mean is efficiency in the formal complexity sense of poly- 
nomial effort (classical or quantum). 

It was shown in Ref. [19] that gapped phases in one dimen- 
sion, in addition to obeying an area law, are actually efficiently 
representable as matrix product states (MPS). The bound on 
the entropy is not very realistic in Ref. [19] (many physical 
systems do exponentially better in terms of scaling with cor- 
relation length) and this bound was improved in Ref. [20]. It 
was also argued in Ref. [18] that even gapless systems may 
be efficiently represented if we are willing to permit a bond 
dimension x for the matrices which grows with system size 
n. Since DMRG can be run in a time that scales polynomially 
with X ix^) and since observables may similarly be measured 
in a time polynomial in x we we see that even if x ~^ poly(n) 
the simulation is still formally efficient [2]. Of course, in prac- 
tice one is always limited in the values of x than can be con- 
sidered, but even two dimensional systems can sometimes be 
handled [21, 22], so the method is quite powerful. 

In higher dimensions it was shown in Ref. [23] that ten- 
sor network states exist with bond dimension x ^ e' 



poly(logn) 



for gapped systems provided a certain bound on the number 
of excited states holds. Of course, this only implies an area 
law up to logarithmic corrections (see also Ref. [24]) where 
as we believe a strict area law is actually valid for these sys- 
tems. Recall that all of our regulated Lorentz invariant models 
in d > 1 obey an area law for the entropy. Depending on how 
easily the network can be contracted, this result may still be 
called efficient. However, in the worst case such a contraction 
is hard to perform, although in such a situation even a con- 
stant X would not be efficient. However, a similar result is not 
available for gapless systems in higher dimensions. It is also 
physically sensible that the leading entropy scaling should not 
change within a phase, and this was recently rigorously shown 
in Ref. [25] for gapped phases. 

Despite the lack of general theorems, we do know of gap- 
less systems or long range correlated systems that have ex- 
act tensor network representations. A simple example is 
provided by a system whose wavefunction is the statistical 
weight of a classical spin model at classical criticality (see 
Ref. [26] for an extensive discussion). If the wavefunction 
is *[cri,....,cr„] = exp(-/3c/i[cri,....,cr„]/2) where ^^ is the 
critical temperature and h is some classical Hamiltonian, then 
correlations of observables diagonal in the classical basis will 
be long-ranged. Furthermore, it follows immediately that the 
system has a tensor network representation and hence that the 
area law is obeyed. 

On the gapped side, we know many examples of topologi- 
cal phases that, in a certain topological limit, have exact tensor 
network representations. These "string-net" phases are topo- 
logical states in d = 2 without protected edge modes. The 
area law is obeyed in all these phases and they have an exact 
and efficiently contractible tensor network representation. Of 
course, we also have models of topologically trivial gapped 
phases that also possess exact tensor network representations, 
e.g. valence bond crystals, at least in a certain limit. It is quite 
reasonable to suppose that for these models a constant bond 
dimension tensor network state is always close to the true 
ground state. A final class of systems worthy of mention is 
the set of non-interacting fermionic topological insulators. As 
shown in Ref. [27] in the non-interacting case (and Ref. [13] 
in the interacting case) the entanglement spectrum of these 
models is quite simple, including the presence of edge modes 
and showing that entanglement is localized near the edge. 



III. LORENTZ INVARIANT ENTANGLEMENT SPECTRA 

We now turn to our main interest, the error in truncating re- 
duced density matrices of spatial subregions in Lorentz invari- 
ant theories. Our tools only permit us to access certain spe- 
cial subregions. For any Lorentz invariant system we have the 
spectrum when _R is a half-space [14, 15] while for any system 
that has conformal symmetry in addition to Lorentz invari- 
ance we also have the spectrum when i? is a d-ball [16, 28]. 
As we indicated in the introduction, these results have had a 
profound impact on our understanding of the physics of black 
holes and quantum gravity. Here we apply the same technol- 
ogy to answer questions of interest in the physics of quantum 



matter 

Consider first the half-space x'^ = t = 0, x^ > 0, and 
x* e R for i = 2, ..., d. As shown in Refs. [13-15, 17] the 
density matrix for this subsystem is given by 



Pxi>o = cxp(-27rJ) 
with J a certain (regulated) boost generator 



J 



x^T" 



(3.1) 



(3.2) 



T°" is a component of the stress tensor and a is a UV cutoff. 
Thus in the notation introduced above we have A'j.i >o = 2ttJ. 
This expression is quite powerful because it gives us the spec- 
trum of the reduced density matrix of a special subsystem in 
terms of a simple geometric problem: the spectrum of the 
field theory on Rindler space ds^ = —p^drf + dp^ + dxf 
(i ~ 2, ...,d). Here the Rindler coordinates i] and p are re- 
lated to the usual coordinates a;° = t and x^ via 



t zL X = ±pe 



±v 



(3.3) 



Rindler space enters precisely because J is the generator of 
time translations in Rindler space, i.e. J generates rj — > i]+Sr], 
so the half-space density matrix is a thermal state of the time 
translation generator in Rindler space at a certain temperature. 

This result has already been used in Ref. [13] to argue 
for a quite general "bulk-edge correspondence" in topologi- 
cal phases that connects entanglement spectra with physical 
edge spectra. This correspondence states that if a physical 
edge has a protected edge mode then the entanglement spec- 
trum will also contain this mode (modulo non-universal de- 
tails) [29, 30]. It is easy to see how Eq. 3.2 implies this result 
since it essentially states that the entanglement Hamiltonian 
looks like the physical Hamiltonian with an edge. More prop- 
erly, these two operators differ by terms localized at the edge 
plus terms in the bulk that don't close the bulk gap. This im- 
plies that any protected edge modes, that is modes that cannot 
be removed by local perturbations at the edge without closing 
the bulk gap, will be preserved in the entanglement Hamil- 
tonian. The ball result discussed below also permitted us to 
give an alternative computation of the topological entangle- 
ment entropy in d = 2 [13]. 

If we interpret Eqs. 3.1 and 3.2 locally, then they give a 
picture of a thermal system with a gap (or energy scale) that 
is growing linearly with x^. Alternatively, we may imagine 
that the temperature is decreasing as 1/x^. The crucial point 
is that, for the half-space, the density matrix is the exponential 
of a local operator in d dimensions (not d—1 dimensions even 
though the entropy scales like the boundary). Furthermore, as 
shown in Ref. [16, 28], the same is also true when i? is a ball 
and the system is conformally invariant. In this case we must 
consider the problem of the field theory not on Rindler space 
but on hyperbolic space H'^ x M (M is time). In both cases we 
know that the entanglement Hamiltonian is a local operator 

Let us use this information to estimate the entropy and re- 
cover the area law for these systems. We focus on the half- 
space region. Consider a conformal field theory (CFT) at tem- 
perature T in d dimensions. The thermal entropy density s 



of this system goes like s ^ T"^. This estimate is always 
valid at sufficiently high temperatures, that is at temperatures 
above all the mass scales in the system. If the system does 
have a gap m then at low temperature the entropy density will 
go like s ^ e~™/^. Now suppose the temperature varies as 
T ^ {x^)^^ as described above. To compute the total entropy 
we integrate the entropy density over the entire half-space 



S 



dx<f-^XA^s{T{x)). 



(3.4) 



This integral must be cutoff at a; ^ a and we assume the trans- 
verse directions x± have size i''^^. The resulting entropy is 



S ' 



d-l 



(3.5) 



which satisfies the area law even though we considered a gap- 
less system (a CFT). If the system has a gap then the x integral 
will be cutoff when x ^ 1/m so that only regions within l/m 
of the boundary at a; = a contribute to the entropy. 

We can do even better when d ^ 1. In terms of the central 
charge c of the CFTi+i the entropy density is 



s(T) = ^T. 



(3.6) 



The temperature is given by T{x) = 1/{2ttx) so that we find 
a total entropy of 



S 



dxsiTix))^-\ogi^/a) 



(3.7) 



where ^ is an infrared cutoff (mass gap). As usual, if the sys- 
tem is truly a CFT then £_ would be replace by a finite interval 
length L. Formally the half space has L = oo so we need an 
additional infrared cutoff. Thus in the one dimensional setting 
we can even get the prefactor of the area law violation correct, 
and yet the entanglement Hamiltonian is still a local operator 
Similar estimates are also possible for the d-ball using the 
connection to hyperbolic space. The metric of hyperbolic 
space may be taken to be 



ds^ 



L\-dt^ 



du" 



smh^udnl_^) (3.8) 



where L is the radius of the ball R and dri^_ 



is the metric 



on the round d — 1 sphere. We must consider the CFT at a 
finite temperature T = 1/{2ttL) on this spacetime to compute 
the entanglement entropy. Note that hyperbolic space is non- 
compact and hence the entropy diverges (as expected). The 
large u cutoff is related to the UV cutoff of the field theory by 
Uc ^ log (L/a). The total volume of hyperbolic space is thus 



V = L'^Qd-i 



dusmh'^-^u^L'^i- 



d-l 



(3.9) 



Since hyperbolic space is homogeneous we have the same en- 
tropy density everywhere, so the total entropy is 



S^sV 



(3.10) 



which again satisfies an area law since s ^ 1/L'^ (recall that 
the temperature is ^ l/L). Furthermore, we can see that when 
d is even there is a universal subleading constant term in S 
while when d is odd there is a universal subleading logarith- 
mic term in S. 

Let us now use these tools to obtain and generalize some 
old results in one dimension. Consider first a critical system, 
i.e. one with conformal invariance. Because of the absence of 
curvature in one dimensional manifolds, the entropy density 
of the CFT on a "hyperbolic line", that is just a finite line 
interval, at temperature T is simply the usual result 



7rc„ 



(3.11) 



Thus the spectrum of the reduced density matrix is simply 
identical to that of a thermal state on a finite interval of length 

V = 2L\n{L/a) at temperature T = {27rL)-'^. Using S = 
sV^ we obtain 



S 



TTC 1 

Y2^ 



(2iln(i/a)) 



-In (L/a) (3.12) 



as required. Moreover, we obtain the complete spectrum in 
terms of the conventional energy spectrum, a result first ob- 
tained in Ref. [31] using the Renyi entropy. Our approach 
has the additional virtue of justifying the weak dependence of 
the non-universal term in S'„ on n (effectively the tempera- 
ture in our langauge). Indeed, we see here that non-universal 
physics enters in the specification of the volume of the hyper- 
bolic space and plays only a minor role. 

The preceeding discussion is sufficient for studying the en- 
tropy, but if we want to ask about individual levels in the en- 
tanglement spectrum, we must specify boundary conditions at 
u = Uc, for example. It can be explicitly checked that differ- 
ent boundary conditions give the same leading scaling of the 
entropy with the central charge, but different boundary condi- 
tions can change the detailed spectrum. A very natural bound- 
ary condition is to demand that there be no energy flux through 
the surface u = Uc (otherwise energy is not conserved). For a 
free scalar field <j>m d ~ 1, the stress tensor is 



The energy current is T"^, so requiring 



T^'lbdy^O 



(3.13) 



(3.14) 



requires dt(t>dx4> = at the boundary. We can satisfy this 
boundary condition by simply setting 9j.0 = 0. These open 
boundary conditions for a scalar field yield a spectrum in 
agreement with the recent numerical calculations of Ref. [32], 
so the field theory methods presented here can provide a 
detailed accounting of the entanglement spectrum of lattice 
models whose low energy limit is the field theory in question. 
We can also obtain interesting results for gapped systems. 
Suppose we now consider one of the above CFTs with a rel- 
evant perturbation that leads at low energies to a finite gap 
in the thermodynamic limit. Because again the entanglement 
Hamiltonian has the form of a local Hamiltonian, it follows 



that states of high entanglement energy can be mapped to 
states of high physical energy. Since states of high entangle- 
ment energy have correspondingly small Schmidt coefficient, 
any attempt to bound the decay of the Schmidt coefficients 
should focus on such high entanglement energy states. The 
difference from the CFT case is that we should take the di- 
mensionless volume of space to be 2 In (^/a) with ^ the cor- 
relation length (recall that L sets the units of all quantities). 
Now the high energy density of states in a CFT in volume V 
at energy E (also made dimensionless) is 



Vl{E) = cxp 



/27rc, 



-VE 



(3.15) 



and the partition function at dimensionless temperature T 

l/(27r)is 



-2ttE 



Z = ^Q.{E)e-'^'' 



(3.16) 



E 



This partition function is just the trace of e~^ and hence is 
proportional the sum of all the Schmidt coefficients. We can 
obtain a sense of how fast the Schmidt coefficients at large E 
decay by arranging the Schmidt coefficients in non-increasing 
order Ai, A2, ... and bounding the function 



f{N)= Y. A, 

i=N+l 

This is nothing but 



-2-kE 



E>E* 



with E* defined by 



N=Y. n{E). 



(3.17) 



(3.18) 



(3.19) 



E<E* 



Since n{E) grows as the exponential of \/E, the relationship 
between A'^ and E* is approximated by 



N = n{E*) 



E 



. 3(ln7V)2 



2ttcV 



(3.20) 



(3.21) 



At large E* or large N we then approximate Z> {E* (N)) by 

3(lniV)2 



/(7V) = Z>(i?*(7V))~cxp 



2clii(C/a) 



(3.22) 



This is a result of relevance for the performance of DMRG 
first obtained for some integrable models in Ref. [33]. Here 
we obtain it much more easily without resorting to integrabil- 
ity and for any Lorentz invariant gapped state. 



IV. GENERAL ARGUMENTS 

We have gone to some effort to establish the intuition that 
these entanglement Hamiltonians may be treated locally. Ul- 
timately our argument is a simple extension of this observa- 
tion. In a nutshell, locality of the entanglement Hamiltonian 
implies that a conventional thermodynamic (large \R\) limit 
will guarantee a concentration of measure around the e"^ most 
probable states. Looking, for example, at the entanglement 
Hamiltonian for the half-space 



K = 2tt 



x'T° 



(4.1) 



we see that the crucial point is that the linearly rising energy 
scale effectively cuts off the entropy in the x^ direction even 
for gapless systems. For a system with a gap our claim is 
intuitively quite reasonable, but below we will show explicitly 
that it is valid even for certain gapless systems. In the case of 
CFTs, the situation is particularly clear when i? is a disk and 
we can use the map to the hyperbolic setting. There we have 
a completely conventional thermal system, albeit in a curved 
space, with a uniform temperature and entropy density. 

Here is yet another perspective. Instead of asking about 
the full spectrum, suppose we refine this information into the 
Renyi entropy 



Sc.{R) = 



1 



1 



■log(tr(p^)). 



(4.2) 



When i? is a d-ball we know that pn has the same spectrum 
as 



Z-^exp{~2TrLHH'i) 



(4.3) 



where Hh^ is the Hamiltonian on hyperbolic space. To com- 
pute the Renyi entropy we consider the more general state 



Zq, ^ cxp [~2-naL'Hijd) 



(4.4) 



with Zq = tr(exp (— 27rQ;LH/fci)). The limit a -^ 00 projects 
onto the ground state ofHfjd while the limit a ^- is a high 
temperature limit. Now at very high temperatures a CFT will 
not detect any small local curvature since high temperature 
corresponds to short distances. Thus we can use the usual flat 
space density of states for the CFT to estimate the free energy 
density. In terms of aFa — — log {Za) we have 



Dq 



aFi — aFa 



1 



a 



(4.5) 



Now we estimate Fa by considering the asymptotic number 
of states n{E) of a CFT in volume V at energy E: 



n{E) ^cxp(k{E'^vy/'-''+^A. 



The partition function at high temperature is thus 



(4.6) 



(4.7) 



which is maximized when 



where L is the Hnear size of R. Thus we must at least have 



1 
T 



dk (V 



1 \E 



E = V 



kdT 



i/(d+i) 



d+l 



(4.8) 



(4.9) 



The free energy is therefore 

F/T = -\ogZ = -V [k 



kdrV If kdT^'^^^'^ 



d+1 



T \d+l 



-V- 



kd 



d+l 



d \d+l^ 
Setting T ~ {2nLay^ we have for F„ 



(4.10) 



aFa = V 



1 



kd 



d{2TTLaf \d+l 



d+l 



(4.11) 



This result is valid at small a since this corresponds to high 
"entanglement temperature". 

Since the conventional flat space thermal entropy is 



__dF 



(4.12) 



we have 



Sth^-{d+l){F/T). (4.13) 

Now in the a — >■ limit we have \aFi \ <C jaFo, |, so 



On 



-aFr, 



(4.14) 



which is indeed positive and is precisely related to the conven- 
tional thermal entropy at a particular temperature and volume. 
In particular, we see that the Renyi entropy goes like (see Eq. 
3.9) 



Sa ^ V{Lc 



\-d 



L a 



(4.15) 



which still satisfies a boundary law. This result is already quite 
interesting since we have shown that the Renyi entropy of a 
ball, as a — > 0, is controlled by the high temperature thermal 
entropy in any CFT (see Ref. [34] for more results for free 
CFTs and Ref. [35] for results about the interacting 0{N) 
model). There are a few subtleties in this result related to 
boundary conditions and topological field theories which we 
discuss in Appendix C. 

Note that the limit a — > is ultimately singular since it 
eventually gives simply the Schmidt rank of the reduced den- 
sity matrix. In any regulated system with a finite number of 
states the growth with a we computed cannot continue for- 
ever At the very least, if there are e'^'^l total states in region 
R (e.g. q = log 2 for qubits) then we must have Sa < q\R\- 
This trivial bound implies that our computation must fail when 



1 



d-l 
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l/d 
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but more generally we expect the CFT result to be valid pro- 
vided the region is large enough and we consider a near one 
so that we are safely probing the long distance physics of the 
regulated field theory. 

Since a < 1 favors the smaller probability states in the tail 
of the distribution, our result for the Renyi entropy can also 
be used to establish that e"^ states is sufficient. Let pM denote 
the density matrix truncated to its M largest eigenvalues and 
let e = Up — Px I It- If we order the eigenvalues of p so that 



Pi 



> 



Pi+i 



then we have e = J2i=M4 



Pi where Sq is the 



Schmidt rank. In Ref. [18] it is shown that the error e satisfies 



log (6) 



< 



1 — a 
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for any < a < 1. Taking 7\/ = (1 — a)e^'''^^ we obtain 

(4.19) 



log(e)<-L^5 



which implies that small error can be obtained by keeping 
roughly e^" states (0 < a < 1). Using our result for any 
CFT that Sa obeys an area law for a not too small, we see 
that again roughly e"^ states suffices. However, if we wish to 
make small error in total system size e ^ l/poly(n) then 5 
must grow with n, however, if we only want small error in 
subsystem size, then we can still respect the A/ ^ e"^ scaling 
independent of ?i. 

To conclude our general arguments, let us discuss the role 
of the regulator in these theories. To have a finite entropy 
we must always regulate the theory, so one might worry that 
the regulator affects our arguments in an important way. We 
argue that this is not the case because the regulator is asso- 
ciated with short distance physics and hence entanglement at 
short scales. The regulator may affect the entanglement of 
nearby degrees of freedom, but its reach is spatially limited, 
so modifying the regulator should only modify short range 
entanglement. This is explicitly visible in our geometric con- 
structions above where the regulator is always associated with 
divergences near the boundary of the region. Formally, chang- 
ing the state inside a large region by adding or removing lo- 
cally entangled degrees of freedom will not change the en- 
tropy of the whole region since these degrees of freedom are 
not entangled with the environment. Such a transformation 
can be modeled by a bounded depth quantum circuit acting on 
the state, but clearly such a circuit cannot qualitatively affect 
the entanglement spectrum e.g. change its entropy scaling. 
This argument is even rigorous for gapped phases up to some 
additional assumptions [25]. However, we emphasize that al- 
though we cannot rigorously prove it, in every calculation we 
can perform all the models we consider have the same scal- 
ing for entanglement through the entire phase whether they 
gapped or gapless. Furthermore, we have a physical renor- 
malization group understanding of this fact. Nevertheless, it 
is interesting to push the rigorous argument as far as we can. 



One way to formalize the idea of local changes to the state 
within a phase is provided by quasi-adiabatic evolution [36]. 
To begin, suppose the phase of interest is gapped and let H{s) 
be a family of Hamiltonians within the phase. We can pro- 
duce a time dependent Hermitian operator D{s) such that the 
ground state |i/'(s)) of H{s) obeys ids\ip{s)) = D{s)\il,'{s)). 
This implies that the ground state of H{s) can be be obtained 
from the ground state at any other point, say s = 0, by evolv- 
ing for a finite time s with the operator D{s). However, the 
crucial part of the story is that D{s) can be a local operator, or 
at least it is local up to (almost) exponentially decaying inter- 
actions. The length scale of the exponential decay is set by the 
inverse gap. What can we say about the entropy generated by 
such an almost local evolution for a finite time (not growing 
with system size)? 

In Ref. [37-39] a bound was estabUshed for the growth of 
entropy of a region A under time evolution with a Hamiltonian 
of the form H = Ha + Hb + Y.i 91^1 <»V^.lf the V' are 
bounded, say 1 1 Vj^ 1 1 < 1 and \\Vg\\ < 1, then we have 



dS{A) 
dt 



<7E 
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with 7 some constant of order one. Note that based on our 
conventions the gj all have units of energy or inverse time. 
Now above we said that the interactions in D{s) were expo- 
nentially localized, so let us suppose that the gj are labeled by 
pairs of point / = {rA,rB) with |g/| < Je^l'"-*^''^!/^. Then 
the sum over / can be written as 
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rAi^A,rB£B 
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but this goes Uke the boundary dA shared between A and B. 

In more detail, if we assume all integrals are effectively cut 
off at scale ^, then every term in D{s) must involve no more 
than ^'^ spins and hence can be decomposed into at most ^ e^ 
product terms of the form Va^Vb- Furthermore, all the terms 
connecting A and B must lie within ^ of the boundary in balls 
of size ^''. Thus we expect at worst a scaling like |9A|e^ ^'^'^^ 
for the entropy bound. Hence time evolution under a Hamil- 
tonian like D{s) for a time of order one will only add at most 
an area law's worth of entropy the ground state. Thus the area 
law is stable within a gapped phase. 

However, it is important to understand that this is not a com- 
pletely rigorous argument since we have imposed a hard cutoff 
at the correlation length. Unfortunately, if we include the ex- 
ponential tail then because of the crude estimate e^ (coming 
from counting the number of ways to decompose an opera- 
tor acting on a region of size ^ into product terms) we cannot 
prove a bound. Now the crude bound we used is almost cer- 
tainly loose, at least for the models we consider which would 
not be expected to have so much entangling power. Indeed, 
the maximum entanglement that could even be generated by 
an operator acting on a region of size ^ is In e^ ^ £,'^ (loga- 
rithm of the dimension of the Hilbert space). If it were true, 
following Ref. [40], that when bounding the rate we could re- 
place the dimension of the Hilbert space on which the entan- 
gling operator acts with its logarithm (as motivated by the fact 



that the maximum entanglement scales logarithmically), then 
a proof would follow. Remarkably, independent of our con- 
siderations here, a claim to a proof of the logarithmic bound 
has appeared [41], so the question may be settled. It would 
then be interesting to extend, if possible, the results of Ref. 
[41] to the Renyi entropy as well (see below). 

Going beyond the entanglement entropy, we can use the re- 
sults of Ref. [37-39] with the entanglement entropy replaced 
by the entanglement Renyi entropy for a < 1 to argue that the 
Renyi entropy is also stable within a phase. Following Ref. 
[37] we maximize 
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with .X e [0, 1] {x and 1 — x are Schmidt coefficients of a 
state in a related two dimensional Hilbert space, see Refs. 
[37, 42]). For 1 > a > 1/2 this maximum exists and pro- 
vides a bound on how fast Sa can grow under time evolution 
with the Hamiltonian Va ® Vb- Combined with the bound 
on the error from Ref. [18] and assuming the same cutoff at 
the correlation length, this argument shows that our truncation 
results are stable within gapped phases. There is one subtlety 
since such a result cannot hold as a — > 0. In that limit the 
Renyi entropy just counts the number of non-zero Schmidt 
coefficients which is a discontinuous quantity. Nevertheless, 
the result is valid for a near one with the breakdown in the 
bound occurring for a < 1/2 since the maximization above 
diverges, i.e. Sa has no finite first time derivative. 

What about gapless phases? Clearly we cannot work be- 
low the gap since then f would be of order the system size 
and nothing could be said. However, there is another natural 
scale in the problem, namely the scale A at which low energy 
physics gives way to high energy physics. For example, in a 
system with emergent photons, A would be the gap to states 
violating the low energy gauge invariance. More generally, A 
is the energy scale beyond which the low energy description 
fails and must be matched onto the high energy description. 
If we now consider quasi-adiabatic evolution with the gap re- 
placed by A then one no longer takes ground states to ground 
states, but we can avoid mixing in states with energy greater 
than A, at least for infinitesimal times. For finite times there 
will be such mixing, but we expect it can be made small so 
that we end up with a state that is not the ground state but 
has nearly the same energy density as the ground state. These 
considerations are similar to those in Ref. [36] in the context 
of emergent gapless gauge fields. The new state will by con- 
struction have no more than an area law's worth of entropy 
compared to the initial state. For example, in a Id CFT the 
leading logarithmic term in the entanglement entropy will not 
be modified by this procedure. In this way we expect that 
gapless phases have the same entropy scaling throughout, al- 
though this conclusion is not completely rigorous. 

Let us mention one final perspective on this issue provided 
by the idea of UV-IR duality (see Ref. [17] for more discus- 
sion). From Eq. 3.2 we see that the states of the entangle- 
ment Hamiltonian of a half-space with large eigenvalue are 
associated not with high physical energy but with long phys- 
ical distances. Even a very low physical energy wavepacket. 



when placed far from the boundary, will have a large expec- 
tation value for the entanglement Hamiltonian because of the 
linearly rising coefficient. Thus it is in some sense the long 
distance physics that is associated with states of high entan- 
glement energy. Such an identification makes physical sense 
since the question of an area law as well as the question of 
approximability should be depend on low energy modes, not 
short distance physics. This is again explicitly visible above 
where the universal parts of the entanglement entropy e.g. the 
log L piece in d = 1, are associated with long distances in 
Rindler space. Since it is the effects of high entanglement 
energy modes that we must control to show that a low error 
approximation exists, it is thus quite reasonable to make such 
calculations within the low energy theory. 



V. ENTANGLEMENT SPECTRUM OF A FREE 
RELATIVISTIC SCALAR FIELD 

We now make our arguments totally concrete by working 
out in detail the case of a half-space in the theory of a free 
scalar field. We first record the spectrum of the entanglement 
Hamiltonian and then give a standard coding argument to ad- 
dress the error 

Following our discussion above, consider a d + 1 dimen- 
sional scalar field in Rindler space ds^ = —p^dif' + dp^ + 
dxf {i = 2, ...,d) (see Ref. [17] for more discussion). The 
scalar wave equation is 



Del) - m^cj) = —^du 
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(5.1) 



Using \J—g = p we find 



1 



2^1^ 



p 



Setting (j> = (t>u{p)e'^-''^ 
ential equation 



0. (5.2) 
"^'' we obtain the ordinary differ- 



^'^4>ui + pdp{pdp4>uj) - p^ik^ + m^ 



0. 



(5.3) 



If we now change coordinates to the "tortoise" coordinate u 
log {p/a) then we obtain the simple equation 



dicp^ + uj'4>u, ~ e""a"(fc" + m")(/.„ = 



(5.4) 



We impose two boundary conditions on this equations corre- 
sponding to the IR u ^ oo and UV u — > (p -> a) limits. 
The IR boundary condition is 4)^{u — > oo) — > 0, and the UV 
boundary condition is <p^{u = 0) = 0. 

We may get a rough idea of the spectrum by replacing these 
boundary conditions with box-like conditions: (t>ui{Q) = 
and 0^(— log (/ca)) = 0. The latter conditions roughly de- 
termines when the "potential" term becomes large and forces 
(f)^ to zero. For each k we then have the standard box "energy 
levels" 



A slightly more refined estimate is provided by the semiclas- 
sical estimate 



Trn 
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log 
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(5.6) 



which reproduces our crude estimate for n <C ka and after- 
wards gives a slightly sublinear growth of w with n. We will 
use the linear spectrum for the remainder of our calculations 
since the results are insensitive to the subleading terms. Al- 
though we will not need them, the exact solutions are simply 
Bessel functions. 

As a simple check let us compute the entropy of this dis- 
tribution for d ~ 1 spatial dimensions (so that k = 0). For a 
boson of frequency oj at temperature T the occupation number 
probabilities are 



p{n) = (1 - e-'3")e-'^'^" 



(5.7) 



and the entropy is 



5 = ^ -(1 - e-'5")e-^'"" (log (1 - e-^'^) - ^ton) . 

(5.8) 
It is convenient to compute this entropy from the heat capacity, 
so we study the average energy of such a collection of energy 
levels. At temperature T (later we set T = {2tt)^^ as in Eq. 
3.1) the average energy is 
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Assuming ma ^ 1 we may approximate the sum by an inte- 
gral as 



(^> 
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The integral may easily be computed to yield 

{J) = ^""^ ^ , 
SO the heat capacity and also the entropy are 
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Using C = TdS/dT we see that in this case C = S, and 
putting T = (27r)^^ we obtain 



S = 



In(^) 

V ma / 
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(5.13) 



which is precisely the entanglement entropy of a nearly mass- 
less scalar field with long correlation length m~^. 

We saw above that the massless scalar field in d + 1 dimen- 
sional Rindler space {d > 1) has a relatively simple spectrum 
that may be approximated as 



(- log ^/ika^^'+Jma)^)^ ' 
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where k is the momentum in the d — 1 transverse directions 
and r = 1, ... labels the Rindler modes. Recall that a is a UV 
cutoff of some sort (a length) and we compactify the trans- 
verse directions in Rindler space into a torus of dimensions 
L X ... X L (L is the subsystem size). The probability to ob- 
serve the occupation numbers Ukr at the Rindler temperature 
27r (see above) is simply 



Pi{nkr}) = Z 



"^exp -27r^Wfernfer (5.15) 



with 
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Note that p factorizes so that different nkr are uncorrected. 

From here on we give a standard coding argument adapted 
to our situation. To that end, let us consider the variable 
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where S is the entropy of p. We see immediately that (y) = 
since (— logp) ~ S. We now compute the variance of y as 
follows. Since the Ukr are uncorrected we have 
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This average may be immediately computed to yield 
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We want an upper bound on the variance of y so using 
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(vahd for a positive) we have 
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with ujk = ^kr/r independent of r. 
Doing the sum over r finally gives 
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We now convert the sum over k to an integral with the upper 
limit |fc| -- 1/a (the UV cutoff) to find 
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where Vld-2 is the volume of the d — 2 sphere. Let us now 
analyze the limits of this expression. First, since ujk -^ oo 
when ka -^ 1 and since the large w limit of integrand goes 
like e~'^^'^ we see that the upper limit of the integral strongly 
converges to zero. Second, the integrand goes like lo^^ at 
small Lu which occurs when fc — > 0. Near the fc = limit we 
must thus consider an integral of the form 

f dkk"^-^ {-log kaf (5.24) 

Jo 

which still converges despite the weak logarithmic singularity. 
Thus at the end of the day the whole integral is bounded by 
some constant C" independent of L and a: 



(,^).C'g) 
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Using this result along with the observation that higher mo- 
ments of y fall off even faster with L permits us to approxi- 
mate the probability distribution for y as 



piy) = 



exp 



2(r^ 
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If we now want the probability that |y| > (5 to be less than e 
we see that asymptotically we need 



6^ . (I 



Thus we conclude that 
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to have small error e. Note that if the variance is paramet- 
rically smaller than our estimate then this condition will be 
achieved even more quickly as a function of L. 

How many states have \y\ < SI The probability of all these 
states is greater than (1 — e) by assumption. Hence we have 



1 > Y, P(Kr}) > 1 - e, 

\y\<s 

but we also know from the definition of y that 

^-S+{L/af-'S > p > ^-S-{L/a)-'''S^ 
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Hence the number N{e, S) of states with \y\ < 6 is bounded 
by 



1 > J2 Pd^kr}) > J2 
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and 



l-e< Ep«"^-»^ Y, ,-S+iL/ar-^S 
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(5.32) 



Combining these inequalities we have 

^s+iL/ar-^s >N>{1- e)e^-(^/'')'"^ (5.33) 

What this finally tells us is that for sufficiently large L we can 
always obtain a good approximation to the density matrix of 
the massless scalar in a certain region of linear size L provided 
we keep roughly e'^ states. 

We have not been very careful about the zero mode in our 
calculation above, indeed, its presence depends on boundary 
conditions anyway. We can always remove this zero mode by 
adding a small regulator mass, and since the contribution of 
zero mode to the variance of y is logarithmic, a very small 
mass (exponentially small in subsystem size) suffices to en- 
sure that the zero mode can be regulated. In most cases such 
a small mass will have no effect on the physics, although one 
place where the zero mode does play an interesting role is the 
case of symmetry broken phases (see Ref. [43]), but this is 
beyond our scope here. 

VI. DISCUSSION 

We have shown that for a wide variety of systems obeying 
the area law, it suffices to keep e'^ states in the reduced density 
matrix to approximate said density matrix to constant error in 
trace norm. Such a result is necessary but not sufficient if 
we want to prove that a tensor network with bounded bond 
dimension x can represent the pure ground state of the whole 
system. To demonstrate these results, we have made heavy use 
of geometric technology available for systems which display 
Lorentz and conformal invariance at low energies. Indeed, 
our whole argument is ultimately based on a renormalization 
group picture of the structure of entanglement in local quan- 
tum systems. Our results do invite some speculation about the 
ability of tensor network states to represent ground states . For 
example, to the extent that we are interested only in universal 
quantities, tensor network states and especially the MERA [7] 
(because of its renormalization group structure) should be in 
the same phase even with bounded x- However, depending on 
how small we wish to make the error, e.g. e ~ l/poly(n), 
we will need to make x scale with n (the system size). For 
example, Eq. 4.19 shows that S ^ \ogn will achieve this er- 
ror bound, but this implies that no bounded x will suffice as 
n — J- oo. It would be very interesting to provide further sup- 
port for these speculations, perhaps by looking at the eigen- 
vectors of the entanglement Hamiltonian. 

Our results also shed light on the physical meaning of the 
entanglement entropy, i.e. that it really does control the ef- 
fective Schmidt rank in some cases. It would be interesting 
to consider "one-shot" versions of quantum information pro- 
tocols, e.g. entanglement distillation in light of our results 
here. For example, the max/min-entropies, which reduce in 



our case to various Renyi entropies, have an operational inter- 
pretation [44]. The single-copy entanglement is the amount 
of entanglement deterministically distillable from the single 
copy of a state and is determined by the a ^ oo limit of the 
Renyi entropy. Our results here suggest that we might be able 
to reliably distill more than 6*00 singlets from a quantum state 
of the type considered above. Finally, the mutual information 
also has a nice operational meaning in the asymptotic setting 
where it represents the total amount of classical and quantum 
correlation between two regions. It would be interesting if this 
result could be strengthened to the non-asymptotic setting for 
regulated Lorentz invariant field theories. 

It is also interesting to observe that the argument around 
Eq. 3.4 estimating the entanglement entropy by integrating 
an appropriate entropy density can be generalized, at a heuris- 
tic level, to many other systems. Below we always suppress 
the dimensionful constant that converts temperature to some 
power of length, e.g. the Fermi velocity. Consider, for exam- 
ple, a scale invariant theory with dynamical exponent z ^ 1 
relating energies to momenta uj ^ k^. The thermal entropy 
density in such a theory is s ~ T'^/^ but in keeping with 
the scaling structure we must also take T ~ 1/x^ instead 
of T ^ l/x. Eq. 3.4 then leads to a very similar estimate 
for the entanglement entropy as in the relativistic case. On 
the other hand, consider a Fermi liquid with thermal entropy 
density s ^ kp~^T and with T ^ l/x (since z = 1 ef- 
fectively). Eq. 3.4 now gives an entanglement entropy of 
roughly {kpLY^^ logi in rough agreement with the known 
result (we used L as an IR cutoff). Finally, even non-Fermi 
liquids give the same logarithmic behavior since although the 
entropy density is now s ^ T^/^ we have again T ~ 1/a;^ 
which yields the same logarithmic estimate. Of course, this is 
only a heuristic since, for example, some z ~ 2 critical the- 
ories, like fermions at the bottom of a quadratic band, have a 
factorized ground state and hence no entanglement. 

In conclusion, even though we have not shown that efficient 
tensor network representations exist, our work does support 
the general notion that these states are in the right comer of 
Hilbert space to describe ground states of local Hamiltonians. 
Furthermore, it is important to perform first principles calcu- 
lations of the type described here since they provide comple- 
mentary information to numerical studies of particular lattice 
models. Of course, even if such tensor networks exist, it may 
still be difficult to work with them, for example, quite apart 
from complexity considerations, the (polynomial) scaling of 
MERA computations with x is unfortunately quite fast. Still, 
the future looks quite bright for these tensor network methods. 
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Appendix A: Comparing quantum states 

Here we briefly review the metrics we consider for com- 
paring two quantum states. There are a number of sensible 
ways we might compare two states, some more stringent than 



11 



others. Let us first discuss a metric that focuses on universal 
physics and then turn to our more stringent quantum informa- 
tion criterion. 

Let H{g,n)he a family of local Hamiltonians with ground 
state |V'(.9, n)) and energy E{g, n). The parameter g is inter- 
preted as a coupling constant which may varied to produce 
different quantum phases of matter while the parameter n is 
a measure of the system size, the total number of sites. The 
thermodynamic limit is obtained by sending n to infinity. The 
states \%}]{g,n — ?► cjo)) are said to be in the same phase pro- 
vided physical observables are smooth functions of g. Often 
one focuses on analyticity of the ground state energy as the 
definition of a phase, however, this is not a particularly physi- 
cally useful distinction. 

A gapped phase is one where a gap in the spectrum of 
H{g,n) persists (does not decay with n) as n ^- oo while 
a gapless phase is one where this does not occur Typically, 
the gap of H{g, n) would go to zero as some power of 1/n in 
a gapless phase. In the cases of interest to us, if d is the spa- 
tial dimension, then the gap will close as n^^/'^. Of course, 
these definitions are only really sketches and complications 
do arise. For example, we may want to restrict n to be even 
in a symmetric gapped phase of spin- 1/2 moments so that the 
ground state can be a singlet. Otherwise, the extra spin- 1/2 
moment will lead to a degeneracy that is not accounted for in 
our crude definition above. 

We emphasize, however, that the limit ?i ^ oo is an artifi- 
cial one. Real physical systems always have a finite effective 
size, that is a finite n. Sometimes the effects of finite n are 
quite important e.g. for fluctuations in mesoscopic systems or 
surface effects, and sometimes they matter little as in the bulk 
of gapped phases. 

Let us now consider two different values of g in the same 
phase. What is the overlap between \il]{gi,n)) and \i>{g2,n))l 
In general we would expect ('0(<?i)|V'(52)) ^ e~". In terms 
of the trace norm we expect ||pi — P2||t ^ 2 — e^" and since 2 
is the largest value possible for normalized states we see that, 
statistically, these two states are maximally distinguishable. 
Nevertheless, both states represent the same phase of matter 
What does this mean? 

Here are two examples. If the system in question was a 
two dimensional gapped topological phase then the topologi- 
cal entanglement entropy would be the same in |-0i) and \1lj2) 
[30, 45]. Similarly, if the system in question was a gapless 
phase or critical point and g was an irrelevant perturbation 
then the long distance decay of correlations, that is opera- 
tor scaling dimensions, would be the same whether they were 
measured in 1-01) or \ip2) [46-48]. Of course, we are simply 
reviewing the completely standard notion of universality. Our 
point is merely that two states need not have high overlap or 
be statistically similar in order to described the same universal 
physics. 

That being said, sometimes we want a more stringent notion 
of equivalent states e.g. for quantum computational tasks. A 
more stringent criterion than the universality notion above is 
provided by demanding small error in trace norm between two 
states. For example, consider the two ferromagnetic states 



and 



|^2> = (v/r^^|t)+aU»^" 



(A2) 



with a <C 1. Their overlap is \/T~~c? and hence their dif- 
ference in trace norm is 



||pi-p2||t = 2v/l-(l-a2)». (A3) 

Requiring that this be small (< e) forces us to choose 



l-eV4<(l-a2)' 



or roughly 



2,2 
na < € . 



(A4) 



(A5) 



iV-i) = I t)®" 



(Al) 



In other words, the states on each site must be close to identi- 
cal with an error going like e/\/n. Such a condition leads to 
nearly identical extensive quantities, not merely close results 
for intensive quantities and universal observables. For exam- 
ple, if we measure the total Sz then we find that these two 
states differ in spin by an amount of order na^ ~ e^. 

A constant e (independent of n) already implies that 
the states being compared have non-vanishing overlap even 
in the thermodynamic limit. Indeed, \\pi — p2\\t = 
2^/1 — |(0i|V'2)p = e clearly gives a finite overlap. Such 
a condition should suffice to argue that the two states are in 
the same phase in the thermodynamic limit. For example, 
a change Sg in the coupling constant can generally be ex- 
pected to degrade the overlap like {ip{g)\i^{g + Sg)) ~ e"""^^, 
so to achieve a constant overlap as n — )- 00 we must send 
Sg — > e^/n. Thus we expect in general that finite over- 
lap corresponds to the same phase with regard to universal 
physics. However, we may be interested in an even stricter 
notion of approximation where the error decreases with n e.g. 
e ~ l/poly(n), and in this case we are really setting the two 
states equal to within the maximum resolution achievable in 
principle on a quantum computer 

For the purposes of extracting universal physics the trace 
norm condition is unnecessarily demanding. Thus we may 
hope that different numerical techniques, especially those like 
tensor network states whose efficiency is based on keeping 
a limited number states, have a reasonable chance of access- 
ing universal information in a quantum phase. On the other 
hand, if it were possible that even small trace norm error could 
be achieved using such a tensor network state, then clearly 
our hope would be fulfilled (although there would remain the 
question of efficiently finding such representations). 

Ultimately, our point of view is the following. The trace 
norm condition is very stringent, but if this condition can 
be satisfied with modest resources in a tensor network set- 
ting, then surely tensor networks can reproduce the universal 
physics of local Hamiltonians. In this way we can provide 
non-trivial evidence, evidence that complements existing nu- 
merical studies, in favor of tensor network descriptions for a 
wide range of phases and phase transitions. 
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Appendix B: A pathological example 

To illustrate why, in principle, we need to worry about ex- 
actly how many states need to be kept, we will now discuss a 
very pathological example. In this section we discuss only a 
classical probability distribution and its associated Renyi en- 
tropy. Since the quantum density matrix can always be diag- 
onalized, we may simply think of this probability distribution 
as the diagonal elements of some density matrix in the appro- 
priate basis. 

The many -body system will be bit strings xi ...x„ of length 
n. The probabilities are 



and 



p(0...0)^{l-q) + ^ 



P(e'*e) = — . 



(Bl) 



(B2) 



In simple terms, with probability (1 — g) we take the string 
0...0 while with probability q we take a random n-bit string 
(including potentially 0...0). Clearly we have ^p = 1 since 



q (2" - l)q 
1 - g + — + ^ — = 1- 



(B3) 



This probability distribution can be understood as a thermal 
state for a very strange Hamiltonian. Let Xi be the value of the 
ith bit and consider the Hamiltonian Hbad = ^Yiii^ ^ -^i^- 



The thermal state at inverse temperature /3 is /? 



-PHtad 



/z 



with 



Choosing (3 so that 



2"-l. 



e/3 + 2" - 1 ^ ' 2" 



(B4) 



(B5) 



we obtain the probability distribution above. Note that Hbad 
is very non-local and the inverse temperature satisfies /3 ^ n. 
Now we compute the Renyi entropy Sa which is 

(B6) 
This expression turns out to be quite delicate in the large n 
limit. For a < 1 then ?i ^- oo we have 

5'„ = n.ln(2) + C'(l). (B7) 

For a = 1 then n — > oo we have 

Si=qnln{2). (B8) 

Finally, for a > 1 then n — > oo we find 

So. = 0(1) (B9) 

All this shows that Sa/n is badly discontinuous at large n, 
indeed, Sa for a > 1 is not even extensive. 



Now we evalute the trace norm criterion above. Restricting 
to the largest x — g^^ states we find that 



\\P-Px\ 



(2" - x)q 



q 



(BIO) 



since 2" ^ 2"''. Thus by taking q close to one we obtain 
an arbitrarily bad approximation. On the other hard, consis- 
tent with the theorem in Ref. [18], the Renyi entropy for any 
a < 1 does tell us how many states we need for a good "ap- 
proximation", but this number is simply all the states! And 
again, the Renyi entropy a > 1 is not even extensive. 

We emphasize that this is a highly artificial and contrived 
example. Nevertheless, it indicates that in principle there is 
something to prove. It is especially intriguing to note that 
while the Renyi entropy for a > 1 totally fails to capture the 
right number of states (even qualitatively) in the above model, 
we will show that these entropies are qualitatively similar to 
the entanglement entropy Si and even to Sa<i in our models. 
This observation is important in the context of numerical work 
since since the Renyi entropy is most easily computed when 
a = 2,3,.... 



Appendix C: The thermal problem on hyperbolic space 

In this section we discuss some of the subtleties associated 
with our result relating the thermal entropy density of a CFT to 
the ground state Renyi entropy of a ball. We obtained this re- 
lation using the mapping of the entanglement spectrum of the 
d-ball to the spectrum of a Hamiltonian of the CFT on hyper- 
bolic space. Because hyperbolic space is homogeneous, it is 
very natural to conjecture that the entropy of a thermal system 
on hyperbolic space is simply the volume of hyperbolic space 
times the entropy density. Unfortunately, the volume of hy- 
perbolic space is infinite and requires regularization, and as is 
often the case, this regularization breaks some of the nice sym- 
metry properties of hyperbolic space. In particular, a bound- 
ary is introduced and boundary conditions should be imposed. 
Often these conditions do not affect the universal properties of 
the entropy provided they obey certain reasonable properties, 
e.g. preserve gauge invariance. This is because the size of the 
boundary at the cutoff Wc = In (L/e), using the metric 

ds^ ^ L^{-df + du^ + sinh^ udnl_^), (CI) 



vol{u = Uc) = ^d-lL'^ TT- 



(C2) 



to leading order in L/e. Thus this boundary size is UV di- 
vergent and hence naively contributes only to non-universal 
quantities. 

Nevertheless, the above argument is a little too quick as the 
following considerations show. Consider the special case of a 
CFT which is also a topological field theory. Not only does 
such a system not care about the size of the metric, it doesn't 
care about the metric at all. Let us think about d = 2 topolog- 
ical states. How can we recover the fact that there is a univer- 
sal topological entropy that doesn't depend on Renyi index (in 
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contrast to the a dependence in Eq. 4.15)? A simple exam- 
ple is provided by C/(l) Chern-Simons theory. This theory is a 
conformal field theory with zero correlation length in the limit 
of infinite gauge field mass, nevertheless, the theory possesses 
gapless edge states on any manifold with boundary no matter 
how large the bulk gauge field mass is. We already indicated 
above that hyperbolic space, when regulated, has a boundary. 
Thus although the bulk is inert, the entanglement spectrum 
(spectrum on hyperbolic space) always has edge states. It is 
these edge states which lead to topological entanglement en- 
tropy. 

In fact, we may explicitly show this as follows [30]. The 
partition function of such 1 + 1 dimensional edge states with 
Hamiltonian H and temperature /?~^ is given by 



Me-'"), 



(C3) 



where a denotes the topological sector (boundary condition); 
we are interested in the identity sector a = 1. A modular 
transformation r — >■ — l/r (r^:, Tj, are the two sides of the 
space-imaginary-time torus) exchanges space and time and 
acts on Z„ as 



Zair) = S'Jbhl/r). 



(C4) 



We already saw above that the temperature is fixed at /3 = 
2ttL, but the size of the boundary is large ~ L{L/e). Such 



a modular transformation is useful because all the modularly 
transformed terms Z have length ^ L and inverse temper- 
ature L'^/e and hence are exponentially close to the ground 
state. Thus we need only keep track of the Casimir energy of 
the identity sector plus the modular S-matrix to compute the 
partition function of interest. At the end of the day the univer- 
sal properties are controlled by the modular 5'-matrix S^^ and 
in particular we find an entropy of In 5";^ = — In V in terms of 
the total quantum dimension V. This was the line of argument 
followed in Ref. [13] to establish the existence of edge states 
in the entanglement spectrum. 

Finally, we see that the other important property of the 
boundary contributions is that they are expected to scale as a 
lowerpowerof "entanglement temperature," e.g. 1/a. Hence 
in the limit of a — > 0, where we claimed that the Renyi en- 
tropy was controlled by the conventional thermal entropy den- 
sity, the boundary contributions are parametrically supressed 
and our bulk argument is valid. More generally, there will be 
subleading in 1/a corrections to this statement; these correc- 
tions are amply demonstrated by the case of topological theo- 
ries which properly have zero thermal entropy (in the limit of 
infinite mass as required to treat them like CFTs) but do have 
a small a independent Renyi entropy. Thus our statement ap- 
plies to the 1/a'' term in the Renyi entropy which is always 
controlled by the conventional thermal entropy even when this 
entropy is zero. 
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